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ABSTRACT 


Normal  mode  theory  is  applied  to  undamped  linear  elastic 
structures  represented  as  lun^>ed  parameter  systems  undergoing 
translational  motion  in  three  directions.  The  derived  equations 
are  primarily  concerned  with  the  re^mnse  of  such  structures 
subject  to  applied  forces  and  base  moticms  and  the  inertia  forces 
required  to  calculate  stress  in  each  mode  of  vibration.  Additional 
relationships  are  presented  for  special  types  of  loading  and  for 
the  effective  mass  acting  in  a  given  mode  due  to  base  motion. 
Similar  equations  are  summarized  in  a^i  appendix  for  structures 
with  six  directions  of  motion,  namely,  three  translational  direc¬ 
tions  and  three  rotational  directions  subject  to  prescribed 
assumptions. 


PROBLEM  STATUS 

This  is  an  interim  report  on  one  phase  of  the  problem;  work 
is  continuing  on  this  and  other  phases. 


AUTHORIZATION 

NRL  Problem  F02-18 
Project  SF  013-10-01-2960  (BuShips) 


Manuscript  submitted  September  Z,  1964. 


SYMBOLS 


A  dot  over  a  symbol  in  the  text  indicates  differentiation  with  respect  to  time. 
d/(  t )  Duhamel  integral  for  base  motion,  direction  r 
F^(t)  applied  force  acting  on  m.,  direction  r 

inertial  force  plus  applied  force  at  in  mode  a,  direction  r 
1 7  impulse  applied  at  m.,  direction  r 
M  total  nuss  of  a  structure 
m/  apparent  mass  in  mode  a,  direction  r 

participation  factor  in  mode  a,  direction  r 
o[^  inertia  force  acting  on  in  the  mode  a,  direction  r 
velocity  step,  direction  r 
x/  absolute  displacement  of  mj ,  direction  r 
xT  relative  displacement  between  and  the  base,  direction  r 

x[^  relative  di^lacement  between  m.  and  the  base,  for  mode  a,  direction  r 
x[^  normal  mode  shape  for  mode  a,  direction  r 
''  base  motion,  direction  r 
itli  mass 

t )  time  function  for  displacement 
t  time 

influence  coefficient 

*  j 

natural  frequency  of  mode  a  for  an  undamped  multi-degree-of-freedom  system. 
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NORMAL  MODE  THEORY  FOR  THREE-DIRECTIONAL  MOTION 


INTRODUCTION 

In  recent  years  normal  mode  theory  has  become  more  widely  used  and  accepted  as 
a  tool  for  structural  design  and  analysis.  While  the  theory  has  been  presented  for  uni¬ 
directional  motion  by  earlier  works,  including  NRL  reports  (1-3),  it  was  felt  necessary 
to  extend  the  theory  to  structures  undergoing  translational  motion  in  three  dimensions. 
Wliile  no  claim  is  made  to  originality  of  the  essential  contents  of  the  report,  many  steps 
are  included  which  often  are  not  published  in  works  dealing  with  the  subject. 

The  primary  concern  of  thiS  report  is  to  find  the  motions  and  inertia  forces  for  cal¬ 
culating  stresses  of  undamped  linear  elastic  structures  which  are  idealized  as  lumped 
parameter  systems  capable  of  undergoing  translational  motion  in  three  directions.  The 
derivation  of  the  equations  is  deliberately  limited  in  the  use  of  mathematical  methods  to 
those  which  are  no  more  complex  than  necessary.  While  this  report  is  self-contained, 
Ref.  3  is  especially  recommended  as  reading  material  to  precede  this  report. 


NOTATION 

2*(vESTiCAL) 

Figure  1  shows  the  orientation  of 
mass  with  relation  to  the  possible 
motions  of  the  base  of  the  structure.  The 
X  axes,  which  describe  the  absolute  mo¬ 
tion  of  a  mass,  are  parallel  to  the  z 
axes,  associated  with  base  motions.  Note 
that  the  origin  of  the  z  axes  is  not  located 
at  a  particular  point;  hence  no  loss  of 
generality  is  made  if  the  1,  2,  and  3  di¬ 
rections  refer  respectively  to  the  longi¬ 
tudinal,  vertical,  and  athwartships  direc  - 
tions  of  a  ship.  It  is  assumed  that  rr.  is 
a  point  mass  with  no  rotational  inertia 
and  that  there  are  n  mass  points  repre¬ 
senting  the  structure. 

Representation  of  Displacements  and  Forces 

The  general  rule  of  notation  for  displacements  and  forces  is  as  follows: 

A  subscript  refers  to  the  mass  point  and  the  superscript  refers  to  the  direction  of 
the  displacement  or  force.  For  example,  represents  the  abrsolute  displacement  of  n 
in  the  rth  direction;  F,'  represents  the  force  applied  at  m.  in  the  s  direction.  Note  that 
1 .  i  1.2.  .11  while  r,  s  1.  2.  3- 

Raising  a  Quantity  to  a  Power 

The  general  rule  of  notation  for  raising  a  quantity  to  a  power  is  as  follows:  Always 
place  the  quantity  inside  brackets  before  raising  to  the  power.  For  example,  to  square 
the  rule  requires  ( 


Z*(AThwaRTSMiPSI 

Fi^.  1  -  Orientation  of  mass  with 
relation  to  the  orientation  of  a  ship 
for  the  case  of  three-directional 
motion 
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Summations 

Unless  otherwise  indicated,  all  summations  on  i ,  j ,  and  k  are  taken  from  1  to  n. 
For  example, 


1=1- 
i  •  •  1 

All  summations  on  a  are  taken  from  1  to  3n.  For  example, 

In 

L  -  L  ■ 

m  sal 


All  summations  on  p,  r,  and  s  are  taken  from  1  to  3.  For  ex:imple, 


L  -  L- 


All  double  summations  are  aU}reviated  in  the  following  form: 


L  -  LL 

i .  f  i  r 


It 


L 

i  “  1 


1 


L  • 


Influence  Coefficients 

The  influence  coefficient  reads  as  follows:  The  deflection  at  i  in  the  rth  direc¬ 
tion  due  to  a  unit  force  at  j  in  the  s  direction.  Thus,  if  a  static  force  is  applied  to  a 
linear  elastic  structure  which  is  attached  to  an  immovable  base,  the  deflection  ^e  to 
distortion  of  any  point  j  on  the  structure  in  direction  r  is  given  by  the  relationsh^) 


rt  t 

L  h,.F,  , 


For  applied  forces  at  each  mass  point  of  the  structure,  this  becomes 


(1) 


Appendix  A  shows  the  influence  coefficient  written  out  in  the  form  to  be  used  for 
finding  natural  frequencies  and  normal  mode  shapes  by  the  iteration  method  (4). 

For  linear  elastic  structures.  Maxwell's  law  of  reciprocal  deflections  (5)  holds, 
namely,  8”  =  5*[. 


FREE  VIBRATIONS 
Normal  Modes 

Assume  that  a  weightless  structure  attached  to  a  fixed  base  is  carrying  a  set  of  n 
concentrated  masses  which  are  attached  at  the  n  points  i.  Consider  its  free  vibrations, 
that  is,  the  possible  motions  in  the  absence  of  external  forces.  This  is  done  by  D'Alem¬ 
bert's  princ4>le,  which  states  that  a  system  in  motion  can  be  cmisidered  to  be  in 
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equilibrium  at  any  instant  if  sqiprt^rlate  inertia  forces  -Bji'  are  applied  to  tiie  system. 
For  the  case  of  the  freely  vibrating  structure,  8itxq>ly  apply  these  inertia  forces  so  as  to 
view  the  structure  as  being  in  a  eiate  of  equilibrium.  The  set  of  forces  on  the  structure 
is  now  treated  as  a  problem  of  statics. 

Recall  that  for  an  elastically  distorted  structure  in  equilibrium 


r  .  ft  • 

X.  =  L  . 


(1) 


For  free  vibrations  the  only  forces  <»i  the  structure  are  the  inertia  forces,  so 

(2) 

This  is  a  set  of  differential  equations  with  constant  coefficients  j^rejts'ng  the  Xj 's 
in  terms  of  the  X  's.  Since  there  is  no  base  motion,  =  X.'  and  f  ■'  =  X^.  Equation 
(2)  is  rewrite 


L  Sji  i«i  Xi 
i  •  • 


*“•  •  • 

s.j  «.  X. 


(S) 


To  obtain  a  solution  try  =  Xj'  sin  (&>t  >  ^),  which  is  usually  dcme  for  the  single- 
degree-of-freedom  system.  Then 


(«)*  £  5 

i .  • 


rt 

ii 


(4) 


Equaticm  (4)  consists  of  three  sets  of  n  algebraic  equations  vdiich  are  written  out  over 
the  range  cm  r  as  follows: 


Xj'  =  (a,)»  I  {i]l  m,x;  +  h]l  m.  Xj*  +  m.  X.') 


L  («”  m.  Xj*  +  S”  n,.  x/  +  S”  m.  X,’) 


X.^  =  L  (s|!  mjX/  +  Sj?  m.  Xj*  +  s”  mj  X.’) 


These  equaticms  can  be  further  written  out  as  3n  algebraic  equations,  i^pendix  A  shows 
these  equations  written  in  matrix  form. 

If  a  solution  is  to  exist  other  than  the  trivial  one  where  all  the  x/'s  equal  zero 
(static  equilibrium  case),  it  occurs  only  for  those  values  of  ui  which  make  the  determi¬ 
nant  of  the  coefficients  of  the  x/'s  equal  to  zero  (6).  This  leads  to  an  algebraic  equaticm 
of  degree  3n  in  usually  called  the  frequency  equation.  Since  undamped  structures 
are  considered,  these  roots  are  real  and  positive  (6).  These  frequencies  are  called  the 
fixed  base  natural  frequencies  of  the  system  oscillating  in  the  absence  of  external  forces. 
For  the  systems  where  the  roots  of  (u>)^  are  all  distinct,  the  ratio  of  amplitudes  of  the 
masses  can  be  found  by  the  back  substitution  solution  of  the  set  of  equations,  v4iich  set  is 
defined  by 


K)*  I  mjXi. 


(5) 
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The  X  ■  a  ar(e  called  the  normal  mode  shapes  and  are  defined  by  Eq.  (5)  for  each  mode  a 
in  each  direction  r . 

Those  systems  which  have  a  pair  or  more  of  equal  roots  are  called  degenerate  sys¬ 
tems.  Other  techniques  for  solving  such  a  set  of  equations  treat  them  as  an  eigenvalue- 
eigenvector  problem,  which  is  a  characteristic  value  problem  with  latent  roots.  For  the 
degenerate  systems,  back  substitution  in  Eq.  (5)  does  not  produce  the  set  of  mode  shapes 
Other  techniques  such  as  matrix  deflation  or  special  forms  of  adjoint  matrices  can  be 
used.  It  is  assumed  that  these  mode  shapes  can  be  found  in  order  to  proceed. 


Orthogonality  of  the  Normal  Modes 

To  establish  the  orthogonality  of  the  normal  mode  shapes,  multiply  both  sides  of  Eq. 
(5)  by  and  sum  on  j  and  r.  This  gives 


..r.r  r  v^rs  -s 

2-  fr-  X  .  X.  =  (  L  X  ,  L  ■  X. 

j  ^  )  j  b  j  a  '  a '  j  j  b  j  i  i  i  a 


L-  s  s  r  r 

IT.  X  ^  L  • .  X., 


-  (  '  ..  J  •  /\  Hi  ■  /V  -  t- 

"  .  5  ‘  '  t  r  ‘i  J  i'’ 


rs  sr  .  , 

Since  J  J  :  • . ; .  Also, 


-s  .Vsr  ►r 

by  Eq.  (5).  Equation  (6)  becomes 


2 

..r  .  r  /  'a\  — s  .  s 

L  ",  ^  (-)  L  \b- 


(6) 


Since  i  and  j  are  now  dummy  subscripts  as  well  as  r  and  s , 


There  are  two  possible  cases:  b  =  a,  or  h  l  a.  When  !>  -  a  the  term  in  the  brackets 
becomes  zero  and  the  summation  becomes 


I 


m 

J 


This  is  a  series  of  positive  terms  which  cannot  be  zero.  When  hi  a,  the  term  in  the 
brackets  is  not  zero,  so  that  the  summation  term  must  be  zero.  This  yields  the  orthog¬ 
onality  conditions 


L 


t  0 


(7) 


L 


m . 
J 


0  . 


a  i  b  . 


(8) 
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Note  that  these  orthogonality  relationships  include  a  double  summation,  that  ia,  the 
usual  summation  of  all  mass  points  as  e]q>erienced  in  the  unidirectional  system  and  the 
additional  summation  over  the  three  possible  directions  of  motion. 


Type  of  Normal  Mode  Solution 

The  distortion  of  the  structure  is  completely  described  if  the  set  of  x/'s  is  found. 
Let  the  time  mode  response  at  point  i  be  X'  .  The  total  re^nse  x/  can  be  found  by 
superposition,  that  is,  ‘  * 


At  each  l  in  each  mode  a  there  is  a  relative  amplitude  of  x[^.  There  must  be  a  function 
which  converts  the  x[^  to  That  is,  a  solution  will  be  sought  in  the  form 


r  — r 

Xi.  =  X 


i. 


so  that 


and 


(9) 

(10) 


If  q,  is  found,  the  free  vibration  problem  is  solved. 
Substitution  of  Eqs.  (9)  and  (10)  into  Eq.  (3)  yields 


By  transposition 


I(<i.  E  .x'.q.)=  0, 

and  by  use  of  Eq.  (5)  this  becomes 


0 . 


Multiplication  of  both  sides  by  m.  x|,,  and  summation  over  j  and  r  yields 


0  . 


There  is  only  one  case  when  the  summation  over  j  and  r  is  not  equal  to  zero:  when 
a  =  b .  The  summation  over  the  modes  is  then  reduced  to 


a, 


=  0  . 


(11) 
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This  has  the  solution 


q.(0) 

q,  =  q,(0)  cos  ^ —  5UI  -,t 


(12) 


Returning  to  Eq.  (9), 


COS 


sin 


(13) 


Initial  Conditions 

Assume  that  the  general  initial  ccmditions  at  t  =  o  are  \  '  -  X^(Q)  and  x/  ::  x.'(O). 
Equation  (13)  yields 

x/(0)  ^  l^[.q,(0).  (14) 


Upon  differentiating  Eq.  (13)  and  introducing  the  initial  condition  on  velocity,  ih^s  yields 

X,'(0)  =  L  x[,  q,(0)  (15) 

The  orthogonality  relationship  can  now  be  used  by  multiplying  both  sides  of  Eqs.  (14)  and 
(15)  by  m.  x[(,  and  siunming  on  i  and  r: 

Xi'(O)  =  I  q.(0) 


Therefore, 


L  m.  x[b  xT(0)  =  L  q.(0)  L  m;  X-^,  x] , 


q.(0)  - 


L  m,  x[.  x/(0) 

1  .  r 

2 


(16) 


q.^0) 


(17) 


Substitution  of  Eqs.  (16)  and  (17)  into  Eq.  (13)  yields  the  complete  normal  mode  solution 
for  free  vibrations: 


y-.  x;.  Z '"jX*.x.’ 


i  (0) 


x;,  L^m.x;.x;(0) 
I  n,.(x*.) 


-  sin  Cl'  t  . 

2  * 


(18) 
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RESPONSE  TO  AN  APPLIED  FORCE 

Consider  a  stnicture  which  rests  on  an  immovable  base,  and  sup[>ose  a  force 
applied  to  m,, ,  is  time  dependent  but  independent  of  structural  reaction.  Using  D'Alem¬ 
bert's  principle  and  influence  coefficients,  the  distortion  of  the  structure  is  described  by 
the  3n  equations 


r  ^  ri  **s  ^  T%  % 

X.  .  -  L 


A  solution  of  the  form 


*i'  =  ?  *i.  0. 


is  sou^t.  Substituting  Eq.  (9)  into  Eq.  (19)  yields 


-  ? x:. .  E  s"  f; 


Transposing, 


?(..x'.*«.i  »;n.x:.)=  zcf.- 


Using  Eq.  (5)  this  may  be  written  as 


(19) 


(9) 


(20) 


Consider  e}q)anding  the  e}q>ression  on  the  right  side  of  Eq.  (20)  into  a  series  of  mode 
shapes.  '  ^et 


Multiplying  both  sides  by  m.  and  summing  on  j  and  r  yields 


The  left  side  is  reduced  by  Eq.  (5),  that  is. 


i  .  r 


(21) 


(22) 


The  summation  over  the  modes  in  the  right  side  of  Eq.  (22)  reduces  to  simply  a  by  virtue 
of  orthogonality.  Thus,  Eq.  (22)  is  rewritten 
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Therefore, 


Equaticm  (21)  becomes 


)k 


Fk 


— r  ^  —1  t 

Xi.Ix,.F, 


(23) 


Since  each  component  of  f,,*  is  independent,  that  is,  the  magnitudes  f,,*  ,  f,,^  ,  and  f,^^  are 
separate  and  independent  of  each  other,  an  expansion  in  s  on  each  side  of  the  equation 
leads  to 


Jk 


(24) 


This  defines  the  influence  coefficient  in  terms  of  the  normal  mode  properties  of  the 
structure  for  translational  motion  in  three  directions. 

Substituting  Eq.  (24)  into  Eq.  (20)  yields 


Transposing, 


Now  the  orthogonality  relationship  is  applied.  Thus, 


Therefore, 
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S 


(25) 


Equation  (25)  is  in  the  form  of  the  equation  of  motion  for  a  single-degree-of-freedom 
system,  thus  having  separated  each  normal  mode.  The  particular  solution  is  written  by 
sight,  using  the  Duhamel  integral  form  for  the  single-degree-of-freedom  system.  Thus 


--^——2  {  3(t-T)rrr. 

Jo 


The  desired  solution  is  then 


X. 


r 


■  a 


— if  [- 


V*  ^ 

F, 


(T)l  sin  ^'  ( t  -  T)  rfT 


(26) 


If  more  than  one  force  is  applied  at  different  points  throughout  the  struct\ire  at  the 
same  time,  say  d  points,  superposition  is  used  to  solve  the  problem.  Since  the  deriva¬ 
tion  assumed  ^e  force  to  be  applied  at  m^,  sum  the  d  applied  forces.  In  this  case  the 
particular  solution  is 


-I 


—  r 

X.. 


r'  r  “ 

■I  Is. 


X,3  F,(T) 


sin  V  ( t -  T)  dX 


(27) 


To  find  the  general  solution  add  the  complementary  solution  represented  by  Eq.  (18) 
to  the  particular  solution  of  Eqs.  (26)  or  (27). 


RESPONSE  TO  BASE  MOTION 

Suppose  a  structure  initially  at  rest  is  attached  to  some  base.  Assume  that  this 
undergoes  a  translational  motion  zV  t)  which  is  a  known  time -dependent  function. 

Consider  the  equations  of  an  elastically  distorted  structure: 


X: 


r  s  s 


Using  D'Alembert's  principle,  this  becomes 


r 


-  L 


T  S 

b  ■  m . 

* )  J 


Since  x/  =  x/  -  Z^  Eq.  (28)  is  written 


where  Xj  represents  the  components  of  relative  acceleration.  Let 


(28) 


(29) 
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r 


‘j. 


and  substitute  this  into  Eq.  (20): 


L  s'*.  Lx*  q 

i..  •»  1  ,  J« 


Q.-  I 


ra  -a 

S .  .  m .  2 

» j  j 


(30) 


With  use  of  Eq.  (24),  the  last  term  of  Eq.  (30)  can  be  expanded  in  its  normal  modes,  that 
is. 


i .  • 


Z 


Substituting  this  into  Eq.  (30)  and  rearranging  terr>s  leads  to 


The  orthogonality  conditions  give 


q,  +  q. 


(31) 


Ibis  equation  is  in  the  form  of  the  equation  of  relative  motion  for  a  single-degree- 
of-freedom  system  if  there  is  a  base  motion  and  no  ajpplied  force.  The  particular  solu¬ 
tion  of  Eq.  (31)  is 


q,  = - ^ - 5  f  r  Z  m.  X*,  Z*(T)j  sin  a*  ( t  -  T)  dT  , 

which  gives  the  relative  motion  of  in  the  r  direction  as 

•  Tp  ° 

The  absolute  motion  of  m.  is 


(32) 


,  •'o 


(t  -T)  dT 


(33) 
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Recall  that  the  origin  of  the  Z  axes  is  not  necessarily  located  at  any  particiilar  ref¬ 
erence  point.  This  is  doe  to  the  form  by  which  base  motions  or  inputs  are  usually  pre¬ 
scribed.  The  inputs  may  be  given  in  one,  two,  or  all  three  possible  directions  of  motion. 
For  example,  siq^se  the  base  disturbance  is  prescribed  in  the  longitudinal  direction 
only.  The  motion  of  m.  in  the  vertical  and  athwartships  directions,  each  of  which  is  per¬ 
pendicular  to  the  longitudinal  motirni,  represents  absolute  motion.  This  agrees  with  Eqs. 
(32)  and  (33),  since  each  equation  is  the  same  for  finding  the  absolute  response  in  the 
vertical  and  athwartships  directions,  while  Eq.  (33)  gives  the  absolute  rcq>onse  for  longi¬ 
tudinal  motion. 


SPECIAL  TOPICS 


Impulse 


Consider  the  impulse  applied  at  mass  k.  The  normal  mode  solution  is 


Xi. 


sin  a).t 


where  l,,',  and  are  the  amplitudes  of  impulse  in  the  1,  2,  and  3  directions,  re- 
spectivoly.  Upon  differentiating. 


— r  v'  — •  » 


cos  . 

2  • 


(34) 


Since  the  structure  rests  on  a  base  and  the  masses  were  assumed  to  be  capable  of  inde¬ 
pendent  movement,  the  velocity  of  m.  must  be  zero  at  t  =  o,  so-that 


—  r  v'  — *  » 


0  . 


Suppose  the  impulse  is  spiled  only  in  the  1  direction,  so  that  the  above  equation 
becomes 


=  0  , 


Since  is  not  zero,  this  reduces  to 


0  . 
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Likewise,  if  the  in^nilse  were  applied  in  the  2  and  3  directions  separately,  there  results 


Z 

z 


0 


0  . 


These  three  erqrressions  can  be  r^resented  more  generally  as 


(35) 


Similarly,  the  velocity  of  the  mass  which  is  struck  by  the  impulse  applied  in  the  1 
direction  is  in  the  1  direction  and  zero  in  the  2  and  3  directions  att  =  o.  There¬ 

fore,  with  reference  to  Eq.  (34), 


zz 


*kB  *ka  ^k 


r 


-2  ~1  1 


r  0 


7^  -it*:.) 


0  . 


Similar  equations  can  be  obtained  for  the  impulse  applied  in  the  2  and  3  directions.  The 
resulting  equations  in  general  form  are 


(36) 


(37) 


Sudden  Motion  of  the  Base 

Consider  the  response  of  a  structure  initially  at  rest  to  a  step  change  in  the  base 
velocity.  Let  ,  z„^  and  comprise  the  components  of  this  step  change.  The  normal 
mode  solution  is 
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Xi 


(38) 


Thus,  if  the  step  change  in  base  velocity  occurs  only  in  the  1  direction,  Eq.  (38)  becomes 


•  r 


Xia 


cos 


At  t  0 ,  the  absolute  velocity  of  each  mass  is  zero,  so  that  the  velocity  relative  to  the 
base  is  -7,‘.  Therefore, 


1 


0 


0  . 


Similar  expressions  can  be  obtained  for  a  step  change  in  base  velocity  in  the  2  and  3  di¬ 
rections.  The  general  equations  which  result  are 


(39) 


1. 


(40) 


Define  as  the  participation  factor 


P 


9 


(41) 
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SO  that  Eqs.  (39)  and  (40)  can  be  rewritten  as 


0  . 

r  t  s 

(42) 

1 . 

r  -  s  . 

(43) 

Now  sum  over 


(43),  so  that 


=  3  . 


Equivalent  Forces  for  Base  Motion 

As  a  final  special  topic  let  t ) ,  in  which  the  foice  on  a  mass  is  propor¬ 

tional  to  that  mass.  Assume  that  such  forces  are  applied  to  each  mass  and  C^t )  is  not 
a  function  of  k .  Then  Eq.  (27)  becomes 


f 


V — f'ri 


-"w  Xwa  C  (T> 


j  sin 


-T)  dT. 


(44) 


This  is  precisely  Eq.  (32),  if  c*(T)  =  Z*(T).  Therefore,  the  displacement  response  for 
many  applied  forces  can  be  converted  to  the  relative  displacement  response  due  to  base 
motion  by  substituting  F^’(T)  =  Z®(T)  and  summing  over  all  k. 


INERTIA  FORCES 
Single  Applied  Force 

In  order  to  calculate  for  stress,  it  is  convenient  to  determine  the  inertia  loadings 
that  the  masses  apply  to  the  structure.  It  has  been  shown  that  each  normal  mode  acts  as 
a  single-degree-of-froedom  system  with  certain  characteristics.  If  the  absolute  accel¬ 
eration  of  each  mass  point  m;  is  found,  the  inertia  forces  can  be  added  to  the  structure 
as  a  loading  by  D'Alembert's  principle. 

Consider  the  case  of  an  applied  force  at  with  no  base  motion.  The  equation  is 


..s  s 


V*  — s  s 

F, 


-  q. 


Solving  for  q„. 
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Since 


then 


Z_J  r  /vM' 


(  „)' 


xl. 


(45) 


Rewrite  Eq.  (45)  using  Eqs.  (35),  (36),  and  (37).  Thus,  for  any  mass  but  the  kth  mass 
(where  the  force  is  applied) 


X/  =  X;  =  -  I  (  X[,  q, 

and  for  the  kth  mass 

-  r  ..  r 
Xk  =  X,  =  - 
k 

The  inertia  loadings  are 

q/  =  L  (u),)2  x[^  ,  (  i  ^  k) 

r  r  — r 

Qjj  =  -  Fj^  ^  ^  Xj^  g  g  . 

These  forces  describe  the  inertial  loadings  for  each  mass  point.  At  ,  there  is  an  ap¬ 
plied  force  fJ.  The  sum  of  the  forces  on  is  the  net  applied  force 

+  fJ  =  Z  X^,  q,. 

The  structure  is  therefore  loaded  in  mode  a  by  a  force  system  of  the  form 

"li  x[,  q,  .  for  all  i  (46) 


E  <‘'a)'  X^',  q. 


where 


‘Ib 


sin  t  -  T)  dT  . 


This  force  system  described  by  Eq.  (46)  may  be  used  to  calculate  stresses  in  the  struc¬ 
ture  for  each  mode. 
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Many  Applied  Forces 

Consider  the  case  where  there  are  many  applied  forces  acting  on  a  structure  which 
vary  as  different  functions  of  time.  'Fhe  equation  is 


k=  1  *  K 


so  that 


^  r  ••  r 

Xi.  =  X,. 


d 

y 

1  t 


*— f  —  li  S 

Xi.  L  IXfc.F, 


!  •.>’  *i.  1. 


X; 


-K) 


~  (lx;,F,’  .  Ixl.r. 

rrP  V  \  *  • 


q. 


Upon  summing  over  the  modes  all  terms  in  the  series  e}q)ression  are  zero  except  when 
r  :  s  and  d  :  i ,  according  to  Eqs.  (36)  and  (37).  Therefore, 


i,'  =  ^  X 


-  -  -  z 


^  Xjg  q^  - 


The  inertia  loadings  are 


Oi'  =  -f/  +  Lv-  x'a)'  "-i  x;^  q.  . 
The  net  force  acting  on  each  mass  is 

o/  +  F;  --  L(  -a)'  -"i  x;.  q.  . 


The  structure  is  therefore  loaded  in  mode  a  by  a  force  system  of  the  form 


(  ‘’b)'*  "i  Xj,  Oa  .  for  all  i 


(47) 


where 


L  .,(x'.)  J. 


-t  r  d 

f  I  I? 

J  L  1  X 


Xu  a 


sin  -’  ( t  -  T)  (IT  . 


This  force  system  described  by  £q.  (47)  may  be  used  to  calculate  stresses  in  the  struc¬ 
ture  for  each  mode  in  the  case  of  many  applied  forces. 
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Base  Motion 

Consider  the  equation  of  motion  in  one  direction  for  a  single-degree -of-freedom 
system  subject  only  to  base  motion  in  that  direction: 

X  f  ( -  -  z. 

The  absolute  acceleration  is  therefore 

X  =  X  ^  Z  ^  -  (o)^X. 

In  an  analogous  manner,  the  absolute  acceleration  in  mode  a  for  the  type  of  structure 
under  investigation  is 


and  the  inertia  forces  are 


Qj.  =  mj  Xj.  =  mj  Xj,  q, 


where 


I  4x1) 


J  Itij  Xj,  Z*(T)j  sin  t  -  T)  dT 


1  r  Tv  s-s  ~i 
=  -  —  I  2-  Pg  Z  (T)  sin  -T)  dT 
•  ■'o  L  «  J 


by  Eq.  (41).  Let 


D  -  U! 

a  I 


r  - 


(T)  sin  ui  ( t  -  T)  dT 


(48) 


(49) 


Equation  (48)  becomes 


The  inertia  forces  in  mode  a,  which  are  the  net  effective  forces  for  calculating  stress, 
can  be  rewritten 


F. 


r  .r  ^  X  s 

^is  "  -  "'i  2.  Pa  Da 


(50) 


EFFECTIVE  MASS  WITH  BASE  MOTION 

To  determine  the  effective  mass  present  In  a  given  mode  of  vibration  for  a  structure 
subject  to  base  motion,  consider  the  net  effective  force  in  mode  a  at  m; : 
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r  r  -  r  »  «  s 

F  -  0  -  m  X  E  P  D  . 

Ifc  t  ^  a  H 

% 

The  total  force  acting  in  mode  a  in  the  r  direction  is 

f  V*  — r  »  s 

F  -  -  L  m.  X  L  P  D  . 

a  I  la  ^  a  a 

»  « 

For  the  single-degtree-of-freedom  system  this  becomes 

F  ::  -MD 

where  li  is  the  mass  of  the  structure  and 


(50) 


(51) 


D 


Z(T)  sin 


.<t  -T)  dT. 


Equation  (51)  is  now  eiqKinded  to  give 


f  ^  1  1  r*  — *“  2  2 

F,  --  -  Lm.  X..P.  D.  -  Lmj  X^.P,  D, 


3 

Da 


(52) 


Since  the  components  of  the  base  motion  are  independent  of  each  ether,  Eq.  (52)  gives 
the  mass  acting  in  the  r  direction  due  to  motion  in  the  s  direction  for  mode  a  as 


-  L  HI;  X^ 


1  a  ^a 


(53) 


For  example,  a  base  motion  in  the  longitudinal  direction  causes  mass  to  act  in  a  mode 
in  each  direction  as  follows: 


M 


1 1 


1 


Pa 


mass  in  the  longitudinal  direction; 


M 


21 

a 


mass  in  the  vertical  direction; 


31 


Z— 3  1 

■"i  Xi.  P, 


mass  in  the  athwartships  direction. 

The  latter  two  terms,  namely  and  m,*,  might  be  called  the  cross-mass  terms. 
It  can  be  shown  that 


m'*  -  m”  ,  r  ^  s  .  (54) 

If  the  masses  are  summed  over  the  modes,  there  results 

M  =  ^  ^  Xjg  Pg  =  S  jn .  M 

i  »  i  * 
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^  L  xj,  p,' 

1  A 

Z  T^.  Z  P  ‘ 


using  Lqs.  (42)  and  (^3). 

Similar  statements  can  be  made  for  base  motions  in  the  vertical  and  ath^-artships 
directions.  The  summation  over  the  modes  of  vibration- can  therefore  be  generalized  as 

m""  M  (55) 

m”  --  h'  -  0.  (56) 


Equation  (55)  indicates  that  the  sum  of  all  the  effective  masses  acting  in  the  r  direc¬ 
tion  due  to  motion  in  the  r  direction  for  the  total  number  of  modes  is  equal  to  the  total 
mass  of  the  actual  structure.  Since  m''  is  always  a  positive  quantity,  a  calculation  of 
the  amount  of  mass  remaining  in  the  hi^er  modes  can  be  made  after  the  lower  modes 
have  been  found. 


SUMMARY 

The  essential  relationships  have  been  derived  for  studying  translational  motion  of 
three-directional  lumped  parameter  systems  based  on  normal  mode  theory.  The  ap¬ 
proach  used  to  develop  these  expressions  is  an  extension  of  an  earlier  repdrt  (3)  on  uni¬ 
directional  normal  mode  theory.  Appendix  B  summarizes  the  equations  of  normal  mode 
theory  for  the  case  of  each  mass  having  six  directions  of  motion,  that  is,  three  transla¬ 
tional  directions  and  three  rotational  directions. 
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Appendix  A 


MATRIX  FORM  OF  LUMPED  PARAMETER  SYSTEMS 


The  relationship  between  normal  mode  shapes  and  fixed  base  natural  frequencies  is 


-  r 


(5) 


The  mode  shapes  in  each  of  the  three  directions  of  motion  are  related  as  follows: 


2 


1  3 
)  i 


^  J 


(Al) 


(A2) 


^ — » 
i 


(A3) 


Each  of  these  expressions  has  a  range  of  n,  so  that  there  are  a  total  of  3r  equations 
represented  by  ^s.  (A1)>(A3).  They  can  be  written  in  matrix  form  as 


1 

-X‘}  = 

(5"liml{X‘>  f 

'.m|{X^^  +  i  i  "m;  '  X^} 

(A4) 

1 

tx^}  = 

(ni){X''  + 

im]{X^}  +  [c”][m]-;XJV 

(A5) 

1  J 

1 

{x^f  = 

[«^')[m){X>} 

(b”] 

[m){X^}  +  [m]{X^} 

(A6) 

where,  for  example, 


-1 

X, 


a 


{X*}  r 


20 


NAVAL  lESfAICH  lAtOtATOlY 


21 


Equations  (A4)-(A6)  can  a)so  be  written  in  matrix  form  as 


This  expression  written  out  in  its  entirety  is  as  follows  on  the  next  page. 


Appendix  B 

EQUATIONS  FOR  SIX-DIRECTIONAL  NORMAL  MODE  THEORY 


NOTATION  AND  ASSUMPTIONS 

1.  The  structure  is  attached  to  a  fixed  base  and  is  represented  by  n  lumped  masses, 
each  mass  being  capable  of  translational  motion  along  three  mutually  perpendicular  axes 
and  rotational  motion  about  each  of  these  axes. 

2.  Each  mass  has  dimensions,  so  that  it  has  rotational  inertia. 

3.  Figure  B1  shows  the  orientation  of  the  ship  by  the  primed  coordinates  ^ixed  axes) 
and  the  axes  of  orientation  for  mass  m.  of  the  structure  (moving  axes).  In  addition  to  the 
usual  translational  motions  given  by  the  1,  2,  and  3  axes,  the  4,  5,  and  6  directions  repre¬ 
sent  the  angular  motions  about  each  axis,  rei^>ectively.  Thus,  x/,  X.^,  and  are  the 
components  of  translational  motion  of  mj ,  while  x^^,  X.^,  and  X  *  are  the  components  of 
rotational  motion  of  m. . 


r(Loogitudinol) 


Fig.  B1  -  Reference  axes  for  orientation 
of  a  ship  and  mass  mj  for  the  case  of  six- 
directional  motion 


4.  For  the  purpose  of  developing  the  equations  of  motion  of  a  body  about  a  fixed 
point,  let  6;‘  =  X*,  0.^  =  Xj®,  and  =  x*. 

5.  Let  be  the  moment  of  Inertia  of  m.  about  its  r  axis  and  l['  be  the  product  of 
inertia  of  m .  about  its  r  and  s  axes. 

6.  It  can  be  shown*  that  the  angular  momenta,  h/,  of  a  body  with  respect  to  its  mov¬ 
ing  axes  are 


*S.  Timoshenko  and  D.H.  Young,  "Advanced  Dynamics,"  New  Y ork:McGraw-Hill  Book 
Company,  Inc.,  p.  332,  1948. 
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1  Ill  I2-I  13-3 

Hi  ^  I,  ,  -  li  i  -  li  i 


2  21'1  22-2  23  ’i 

Hi  =1.  i  -  li  i  -  li  . 


3  31  •  I  32-2  33-3 

Hi  =  I.  i  -  li  .  -  li  . 


7.  Select  the  1,  2,  and  3  axes  of  m.  as  the  principal  axes  so  that  the  products  of 
inertia  are  zero  and  Eq.  (Bi)  reduces  to 


1  1 1  •  I 

H.  .  I,  . 


2  2  2  •.  2 
Hi  --  I.  . 


3  33  -3 

H.  =  1;  - 


8.  The  principle  of  angular  momentum  states  that  the  rate  of  change  of  the  angular 
momentum  of  a  body  rotating  about  a  fixed  point  is  equal  to  the  moment  of  all  forces 
acting  on  the  body  with  respect  to  the  same  point.  After  taking  into  account  the  rate  of 
change  of  the  angular  momenta  with  respect  to  the  1,  2,  and  3  axes  and  the  fact  that  the 
1,  2,  and  3  axes  are  also  rotating  about  a  fixed  point,  the  following  equations  result; 


♦  2  *  3 

•  3  2 

=  N. 

H  “ 

H 

dt 

1  1 

1  I 

1 

2 

dH, 

•  3  1 

•  I  3 

i  Hi  - 

i  Hi 

:  N 

dt 

dH^ 

*  1  2 

’  2  1 

H;  - 

.  Hi 

:  N 

dt 

S  I 

where  Nj  is  the  moment  of  all  forces  acting  on  the  body  about  the  r  axis.  Substitute  Eq. 
(B2)  into  (B3)  to  get 


d,!. 

j 

f  33 

2  2  ^ 

,  •  2  ■  3 

1 

+  1 

,  I  . 

I  ■  , 

--  Ni 

dt 

^  I 

»  1 

'  1  > 

•  2 

(  '> 

3  3^ 

1  •  I  •  3 

2 

-  * 

+  1 

I  •  - 

I.  ] 

!■.  r. 

--  N; 

dt 

\  I 

1  1 

dfV^ 

* 

^22 

•  0 

1  •  I  •  2 

3 

+  1 

ll  - 

I  ■ 

1  H. 

r  N; 

dt 

i  1 

>  1 

1  1 

These  equations  are  called  the  Euler  equations  tif  motion. 

10.  For  small  oscillations  assume  that  the  terms  containing  the  product  of  the  o's 
are  small  compared  with  the  other  terms  in  Eq.  (B4)^  so  that 
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111 


N. 


31  1 


I: 


=  N. 


3J  V  3 


=  N; 


Using  the  notation  of  Fig.  B1  for  m. ,  these  equations  become 


4  4  4 


5-5 

I;  X:  =  N; 


6-6 

I:  X:  =  N; 


{B5) 


where,  for  example,  I;*  is  the  moment  of  inertia  of  in  the  4  direction  (or  about  the 
1  axis)  and  is  the  angular  acceleration  of  in  the  4  direction.  Equations  {B5)  rep¬ 
resent  the  relationships  between  the  inertia  torques  and  the  applied  torques  for  equilib¬ 
rium  about  a  fixed  point.  These  inertia  torques  along  with  the  inertia  forces  will  be  used 
for  the  free  vibration  problem  of  the  structure  under  investigation. 

11.  Since  the  axes  of  each  mass  are  the  principal  axes,  the  axes  of  different  masses 
are  therefore  not  necessarily  parallel  to  each  other.  This  requires  a  new  definition  of 
the  influence  coefficient  ‘i*  as  the  motion  (deflection  or  rotation)  of  m-  in  the  r  direc¬ 
tion  of  rtij  due  to  a  unit  load  (force  or  torque)  at  m.  in  the  s  direction  of  m..  For  exam¬ 
ple,  bjj  is  the  deflection  of  m.  in  the  1  directicm  of  m-  due  to  a  unit  force  at  m.  in  the  1 
direction  of  mj.  Note  that  the  1  direction  associated  with  n-j  is  not  necessarily  parallel 
with  the  1  direction  of 

12.  It  is  assumed  that  the  change  in  geometry  of  the  structure  is  small  during  its 
dynamic  response  under  the  action  of  external  forces  and  torques.  This  means  that  the 
influence  coefficients  calculated  for  statical  loads  on  the  structure  with  respect  to  the 
principal  axes  of  each  mass  are  used  to  find  the  dynamic  motions  while  the  axes  are  per¬ 
mitted  to  translate  and  rotate  with  each  mass. 


13.  Define  the  direction  cosine  between  the  r  axis  and  r'  axis  at  mj  as 


r  r  ‘ 

t 


cos^{  r.  r  '  )  . 


Assume  that  during  the  structure's  dynamic  response  the  direction  cosines  remain 
constant. 


NORMAL  MODE  EQUATIONS 

The  equations  of  motion  for  the  free  vibrations  of  the  assumed  structure  are  written 
using  D'Alembert's  principle  and  treating  the  inertia  forces  and  torques  as  the  applied 
loads. 
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(  =  -  E  E  m  X.  -  r  r  6. .  I  X. 

'j  ^  j  1  I  ^  ji  *1 


I  S  ‘  I 


L 

I  t‘  I 


»  rt  %  ‘‘  % 

L  h  Xi  . 


r  M. 


where  Ij  =  mj  for  s  =  1,  2,  and  3. 


(B6) 


Equation  (B6)  is  precisely  the  same  as  Eq.  (3)  except  that  the  range  on  the  direction 
of  motion  includes  six  independent  coordinates  for  each  mass.  The  normal  mode  equa¬ 
tions  r  this  case  are  now  summarized  from  earlier  results  of  this  report,  where  sum¬ 
mations  on  r,  s,  and  p  are  from  1  through  6  unless  otherwise  indicated  and  the  summa¬ 
tion  on  a  is  from  1  to  6n; 


Mode  Shapes  and  Natural  Frequencies 


Xj,  -  I  i-i  l.  Xj,. 


a=  1.  2. 


6n 


(B7) 


Orthogonality  of  Normal  Modes 


1 1; (*:.)’  i  0 
I  Ij  x^.x;,  =  C,  a^  b. 


(B8) 

(B9) 


Response  for  Free  Vibrations 


X,  =  ;  — Lj - — CO.  ..t  *  ;  - 


—  r  ^  s  —1  •  s 

Xi,  ^  1 3  X.,  Xj  (0) 


sin  uj  t.  (BIO) 


Response  to  an  Applied  Force  and  Torque  at 


=  T  — (  f  i  K.  Km  .  i  K,  n; 

^  B  KK)  J.  ^ . 


sin  <»  (  t  -  T)  dT 


Response  to  Many  Applied  Forces  and  Torques 

Assume  that  there  are  d  applied  forces  and  h  applied  torques;  therefore, 


(BID 


r  \  *  X^  r*  r  hi 

/  =/  - — - i  I  ^  X;3F;(T)  f  E  L  sin,Vt-T)dT.(B12) 

L  I.”  x"  )  I  v-l  .M  J 

*  ■  i.p  *  \  *•/ 
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Response  to  Base  Motion 

Consider  the  case  where  translational  base  motion  is  the  prescribed  input  referred 
to  the  axes  of  the  ship,  namely,  the  primed  axes  shown  in  Fig.  Bl.  If  each  mass  of  the 
structure  is  loaded  with  the  set  of  forces  -m^z* ',  the  reqx>nse  of  each  mass  is  the  same 
as  for  the  case  of  many  applied  forces.  However,  this  re^[>onse  is  now  the  relative  mo¬ 
tion  between  each  mass  and  the  base. 

Figure  B2  shows  the  forces  acting  on  oriented  with  respect  to  the  primed  refer¬ 
ence  axes.  To  transform  these  forces  along  the  principal  axes  of  m. ,  use  the  direction 
cosines  between  the  axes  r  and  r'  at  m. .  The  transformed  forces  at  m.  in  the  r  direc¬ 
tion  are 


Fig.  BZ  -  Forces  acting  on  idj  for  finding 
the  response  due  to  base  motion  for  the 
case  of  six-directional  motion 


Use  this  summation  to  represent  the  forces  acting  on  the  structure.  With  reference  to 
Eq.  (B12),  the  relative  motion  between  m.  and  the  base  in  the  r  direction  is  therefore 


f 


sin  t  -  T)  dT  • 


(B14) 


Special  Topic  —  Sudden  Motion  of  the  Base 

Consider  the  response  of  a  structure  initially  at  rest  to  a  step  change  Ui  the  trans¬ 
lational  velocity  of  the  base.  It  follows  from  Eq.  (BI4)  that 
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3  ' 


X'.LZ  £  ..-X: 

It  t-  I  »  '  ^  I  ' 

I  .P  '  • 


—  s 
ka 


COS 


:3i5) 


where  ,  z^  ,  and  z^  represent  the  amplitudes  of  the  base  translational  velocity  in 
the  r,  2',  and  3'  directions,  respectively.  The  six  components  of  velocity  defined  by  Eq. 
(B15)  are  now  treated  as  two  separate  groups,  namely,  the  three  translational  components 
and  the  three  rotational  components. 


The  translational  components  of  velocity  response  are  treated  first.  At  t  =0,  the 
absolute  translational  velocity  of  each  mass  is  zero.  To  find  the  initial  relative  transla¬ 
tional  velocity  of  each  mass  with  respect  to  the  base  motion,  transform  the  base  motion 
into  the  direction  of  the  principal  axes  of  each  mass  as  follows: 


Xj  (0) 


(B16) 


Substitute  Eq.  (B16)  into  the  left  side  of  Eq.  (B15)  for  t  =  o: 


J  ' 


3  ' 


2-1  ^  "'k  x^ 

k  1  1  '=1 ' 


a  1( 


(B17) 


Since  the  components  of  the  base  velocity  are  independent  of  each  other,  it  follows  that 


Let 


r=  1.  2,  3  . 


SO  that  Eq.  (B18)  becomes 


t 


t 

r  r 

i 


I  x;  z  I  p:. 

•  '  k  » »  1 


r  =  1.  2.  3  . 


(B18) 


(B19) 


(B20) 


In  the  case  of  rotational  motion,  the  rotation  of  the  base  is  zero,  so  that  the  initial 
relative  rotational  velocHy  of  each  mass  is  also  zero;  that  is, 


x/{0)  =  X^(0)  =  0  . 


r  ^  4,  S,  6  . 


(B21) 
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Referring  to  Eq.  (B15)  at  t  o, 


Since  the  components  of  the  base  velocity  are  independent  of  each  other,  this  becomes 


J 


*  s  ! 


r  4.  5-  6:  s'  r.  2^  3' 


(B22) 


using  Eq.  (B19). 


STRESS  CALCULATIONS 

To  calculate  stresses  in  the  structure,  apply  the  net  effective  forces  and  torques  at 
each  mass  point  for  each  normal  mode  of  vibration.  With  these  loads  acting  in  each 
mode,  the  stresses  can  be  calculated  for  each  mode  throughout  the  structure,  and  the 
final  stresses  are  obtained  by  superposing  over  the  modes.  This  approach  is  the 
same  as  followed  earlier  in  the  case  of  translational  motion  in  three  directions  only. 
The  net  effect  of  forces  and  torques  for  special  cases  are  summarized  as  follows. 


Single  Applied  Force  and  Torque  at 


(B23) 


where 


q.. 


1 


I 


p 

I 


X.,  F,(T) 


+ 


6  _  1 

Xl,  N;,T)J  s.n 


.,(t  -T)  clT. 


Many  Applied  Forces  (rl)  and  Torques  (h) 


Fin  X,„q„ 


(B24) 


where 


■  I  i  h  6 

^  X^,  F^(T)  *  I,  Z,  N^(T) 

k  l  s  I  k  -  I  4 


Sin  J  t  -T)  ciT  . 


Base  Motion 


(  xin  n, 


(B25) 
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tdiere 


=>.  =  - 


v:r-^-- - :i  f  III  ^  ,  «k  K.  z’  ^T)l  sin  a,,(t- 

/  -  P/kP  \  J  L  k  *=<1  .  -r  J 

fi:  •'o 


T)  dT  ■ 


EFFECTIVE  MASS  WITH  BASE  MOTION 


To  find  bow  much  effective  mass  is  acting  in  each  normal  mode  of  vibration,  first 
consider  the  effective  force  in  naode  a  at  m-  from  Eq.  (B25): 


Fj,  =  «j  X[.  q.  .  r=  1.  2.  3 


where,  using  Eq.  CB19), 


in  vdiich 


'  « '  k  •  1  « '  -  1  • 


D. 


(T) 


sin  t  -  T)  dT  . 


(B26) 


(B27) 


(B28) 


Tbese  forces  are  transformed  to  the  primed  axes  (orientation  of  the  ship's  motion)  as 
follows: 


L< 


Fi 


r'  =  1',  2'.  3' 


(B29) 


Substitute  Eq.  (B26)  and  (B27)  into  Eq.  (B29).  This  gives 

3  3  3' 


r'  ^  .r'r  _r  V’  .»»'  a  »' 

Fi.  =  -  -"i  L  L  L  -tfc  Pfc.  D. 

r=  I 

The  effective  force  acting  in  mode  a  is 

<  =  -  L  L  <■'  E  1 1  <;■  K.  <  •  (“®> 

■  r=l  k  i*!*'*!' 

SO  that  the  effective  mass  acting  in  the  r '  direction  due  to  motion  in  the  s '  direction  is 

“l’’’  -  L  "it  xj.  L  t  p;. . 

i  T-l  k  t=l  " 

The  cross-mass  terms  are  again  symmetrical,  that  is. 


k  •»! 


3  3 
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'M' 


Consider  the  summation  ol  the  effective  masses  over  the  modes: 


(B32) 


-  L<'  =  I*. EC 


i  r»  1 


k  •-  1 


With  reference  to  Eq.  (B20)  this  reduces  to 

,  .  j  ,  , 

r  ■  ,r  r  r* 

“  =  E  “i  I.  ^  i 


(B3S) 


i  I 


It  has  been  shovm*  that 


f  •  1 


=  0  .  r'  i  s'  . 


Since 


<  '  =  t'' 


Eq.  (B33)  becomes 


M*^*  =Z.n>i  =  M.  r'  =  s' 

=  0  .  r'  ^  s'  .  (B34) 

As  in  the  case  of  three-directional  motion,  the  total  mass  acting  over  the  modes  in  the 
r '  direction  due  to  a  base  motion  in  the  r '  direction  equals  the  total  mass  of  the  struc- 
txire,  while  the  summation  of  the  cross-mass  terms  equals  zero. 

The  torques  due  to  translational  base  motion  are  now  treated.  It  is  necessary  to 
introduce  new  notation  for  the  angular  acceleration  and  the  moment  of  inertia.  In  addi¬ 
tion,  matrix  notation  will  be  used  in  transforming  the  torques  from  the  principal  axes  of 
each  mass  to  the  ship's  orientation. 

For  angular  acceleration,  let 


**3 

Xi.  = 


(B35) 


'This  is  necessary  since  direction  cosines  are  used  which  refer  to  the  1,  2,  and  3  axes 
and  the  1',  2',  and  3'  axes.  Likewise,  n/,  and  Nj*  now  represent  the  torques  at 
about  the  1,  2,  and  3  axes,  respectively. 

The  direction  cosines  are  written  in  matrix  form  as  follows: 


♦Herbert  Goldstein,  "Classical  Mechanics,"  Cambridge: Addison- Wesley,  p.  98,  1950. 
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1  1 

2 

13 

}  1 

1  1  , 

^  ■> 

1  “  " 

(B36) 

3  15 

3  2 

3  3, 

When  using  the  direction  cosines  for  transiornulions,  these  transformaiions  arc  called 
ortnogonal  transformations,  and  it  can  be  shov^-n^  that 


(B37) 


and 


(B38) 


7'his  last  equation  relates  the  inverse  niatrLx  to  the  transjxjse  ni.ilrLx  of  the  direction 
cosines. 

The  transformation  of  the  torques  due  to  base  translational  motion  can  now  be  made 
as  follows: 


S' 


s  1.23  N  12  3 


The  torques  referred  to  the  unprimed  a.xes  are  written  in  matrix  form  .as 


where 


1'  .  '  . 


I  ! 

1  (' 


i 

I  0  I 


0  0  I 


Premultiply  both  sides  of  Eq.  (B40)  by  ^ ,  which  gives 


•N'' 


1'  ' 

I  J  I 


i  1  1  I  t 


(B39) 


(B40) 


(B41) 


(B42) 


In  general,  the  equation  for  the  primed  system  relating  the  torque  with  the  inertia 
terms  and  angular  acceleration  is 


fRobert  L.  Halfmaii,  "Dynamics, "  Vol.  1,  Reading,  Massachusetts; Addison- Wesley,  p. 
207,  1962. 
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(B43) 


Comparing  Lerms  in  Eq.  (B42)  and  (B43)  and  using  Eq.  (B38)  and  (B39)  leads  to 


(B44) 


(B45) 


Equations  (B43),  (B44),  and  (B45)  are  now  written  in  the  equivalent  series  form  as 


J  ' 


»  ’  1  ' 


(B46) 


3  3,, 

« — •  c — •  r  u  s  V  u 

--  L  L  'i  ^ 

u  -  1  •  I 


(B47) 


3 

T 
%-  1 


s  s  s 
t  i 


Rewrite  Eq.  (B46)  for  the  mode  a  as 


or 


where 


T*  ^  r*  s'  s' 

-  -  L  I.  .a 


N. 


r  '  s '  s'  s  “S 

L  I.  L  ", 


s-  1 


••  s 
i  a 


-  (  ‘a)"  *1. 


s  1  2.  3 


1  Y"  V"  s's  s  s' 

'  a'  k  s=  I  I  ' 


CB48) 


(B49) 


(B50) 


(B51) 

(B52) 


It  should  be  noted  that  x,'p,  xj^,  and  X- ^  in  Eq.  (B51)  are  actually  x^^,  x^,  and  x* 
respectively,  due  to  the  new  notation  for  the  transformation  on  the  torques. 

Substitute  Eq.  (B51)  and  (B52)  into  Eq.  (B50)  to  get 


3' 


z; 


X: 


,  '^p 


Da 


(B53) 


r  ^  1  p 


w 
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Sum  over  all  the  mass  points  to  get 


1 1  L 


r  =  1  p  ' 


(B54) 


From  Eq.  (B54)  the  effective  moment  of  inertia  in  the  r '  direction  in  mode  .1  due  to  a . 
base  tianslational  motion  in  the  p'  direction  is 


L  L  z;  E 


rp  r 

'  P 

k 


(B55) 


Sum  over  the  modes  to  get 


r  P 


L  E  E  I  c  1 1 


k 


ka 


14*1*  a  kr“I 

It  has  been  shown  that  for  rotational  motion  due  to  a  sudden  translational  motion  of  the 
base  j 


E  E  E  -> 


^  r  rp 


k  a  k 


=  0  . 


(B22) 


k  <  ^  1 


SO  that  /  :  0  . 


For  the  case  where  the  principal  axes  of  each  mass  are  orientated  parallel  to  the 
ship’s  reference  axes,  Eq.  (B55)  reduces  to 


(B56) 


Note  the  similarity  between  this  equation  and  Eq.  (53).  It  has  been  previously  stated  that 
^ile  the  range  on  r  in  Eq.  (B56)  is  1,  2,  and  3,  the  terms  are  in  fact  x^,,  x*,,  and 
X*,,  respectively.  Therefore,  Eqs.  (42)  and  (43),  which  were  applied  to  Eq.  (53)  after 
summing  over  the  modes,  cannot  be  used  with  Eq.  (B56). 
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